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CHAPTER  I 


INTRODUCTION 

In  the  design  and  use  of  large  composite  panels  in  structures  such 
as  aircraft  and  space  vehicles  a  major  concern  is  the  ability  of  the 
panel  to  function  (continue  to  carry  a  substantial  part  of  its  design 
load)  after  being  damaged.  One  technique  that  has  evolved  from  the 
stringer  reinforced  metallic  panel  is  the  buffer  strip  or  hybrid  panel 
shown  in  Figure  1.  Because  of  the  fabrication  methods  used  in  composites 
it  is  possible  to  make  such  a  laminate  by  replacing  specific  fibers, 
usually  parallel  to  the  load  axes,  with  fibers  of  the  appropriate  physical 
and  geometric  properties  necessary  to  arrest  a  crack  that  originates  in 
the  parent  laminate  material.  Since  the  buffer  strips  are  usually  narrow 
and  relatively  far  apart  the  stiffness,  weight,  and  strength  of  the 
undamaged  laminate  is  not  significantly  effected  by  the  replacement. 

Much  experimental  work  has  been  done  to  investigate  this  behavior, 
i.e.  to  determine  the  best  buffer  strip  material,  with  the  studies  of 
Eisenmann  and  Kaminski  [1],  Hess,  Huang  and  Rubin  [2J  ,  Avery  and  Porter  [3], 
Verette  and  Labor  [4],  and  Poe  and  Kennedy  [5]  being  significant  contribu¬ 
tions  in  this  area.  The  same  cannot  be  said  for  published  analytical 
solutions  and  it  is  this  question  that  is  considered  in  the  present  study. 

Some  of  the  first  work  in  modeling  a  uni -directional  composite  con¬ 
taining  broken  fibers  was  presented  by  Hedgepeth  [6]  where  the  case  of 
no  additional  damage  other  than  the  initial  notch  was  considered.  This 
study  was  extended  by  Hedgepeth  and  Van  Dyke  for  the  special  case  of  one 
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Figure  1.  A  Typical  Buffer  Strip  Laminate  Configuration 


broken  fiber  with  matrix  yielding  parallel  to  the  fiber  [7]  and  for  one 
fiber  with  longitudinal  splitting  in  the  matrix  [8].  Goree  and  Gross 

[9]  extended  the  Hedgepeth  solutions  to  include  longitudinal  matrix 
yielding  and  splitting  for  an  arbitrary  number  of  broken  fibers.  The 
results  of  Goree  and  Gross  gave  very  good  agreement  with  experimental 
results  for  brittle  matrix  composites  which  exhibit  large  longitudinal 
matrix  splitting.  For  ductile  matrix  composites  such  as  boron/aluminum, 
which  exhibit  large  yielding  but  very  little  splitting  in  the  matrix, 
this  model  predicted  the  right  trend  but  the  agreement  was  not  very 
good,  especially  for  short  notch  lengths.  Goree,  Dharani  and  Jones 

[10]  attempted  two  modifications  to  the  above  solution.  First,  the 
matrix  was  assumed  to  be  strain-hardening  and  secondly,  a  cover  sheet 
was  included  over  the  main  laminate.  The  results  of  [10]  showed  that 
the  inclusion  of  either  a  strain-hardening  matrix  or  the  addition  of  a 
cover  sheet  did  not  improve  the  agreement  between  the  predicted  and 
experimental  results.  Based  on  the  observed  fact  that  in  addition  to 
longitudinal  yielding  of  the  matrix,  a  certain  amount  of  stable  trans¬ 
verse  extension  of  the  initial  notch  under  increasing  applied  load 
takes  place,  Dharani,  Jones  and  Goree  [11]  then  extended  the  solution  of 
[9]  to  include  transverse  damage  ahead  of  the  initial  notch  in  addition 
to  the  longitudinal  matrix  damage.  The  results  of  [11]  showed  a  very 
significant  improvement  in  the  ability  of  the  model  to  represent  the 
behavior  of  a  ductile  matrix  composite. 

In  all  these  analytical  studies  the  laminate  is  modeled  as  a  two- 
dimensional  region  having  a  single  row  (mono-layer)  of  parallel,  identi¬ 
cal,  equally  spaced  fibers,  separated  by  matrix.  The  damage  is  taken  to 
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consist  of  an  arbitrary  number  of  broken  fibers  such  that  all  breaks 
lie  along  the  x-axiSj  but  they  need  not  form  a  continuous  break  (notch). 
The  fibers  are  assumed  to  be  of  much  higher  strength  and  extensional 
stiffness  than  the  matrix  and  all  the  axial  load  is  assumed  to  be 
carried  by  the  fibers,  with  the  matrix  transfering  load  by  shear 
stresses  as  given  by  the  classical  shear-lag  assumption.  One  very  impor¬ 
tant  feature  of  the  shear-lag  assumption  is  that  is  simplifies  the  equili¬ 
brium  equations  by  removing  the  transverse  displacement  dependence  from 
the  longitudinal  equilibrium  equation.  The  fiber  stress  and  matrix  shear 
stress  can  then  be  determined  without  solving  the  transverse  equation. 

The  methods  of  analysis  developed  and  discussed  in  the  above 
studies  [6-11]  are  extended  in  the  present  work  to  determine  the 
stresses  and  displacements  in  a  hybrid  uni -directional  laminate  having 
an  initial  notch  in  the  vicinity  of  a  single  finite  width  region  of 
different  material  properties  as  shown  in  Figure  2.  This  geometry  is 
an  idealization  of  the  usual  periodic  placement  of  buffer  strips  as 
indicated  in  Figure  1  and  it  is  assumed  that  the  stresses  near  the 
notch  tip  and  the  single  buffer  strip  are  approximately  the  same  as 
those  in  a  wide  panel  with  relatively  narrow  buffer  strips. 

Of  particular  interest  is  the  investigation  of  the  behavior  of  the 
laminate  as  a  function  of  the  relative  ultimate  stress  and  extensional 
modulus  of  the  buffer  strip  fibers,  the  buffer  strip  width  and  initial 
notch  length. 

As  an  initial  step  in  understanding  the  basic  mechanism  of  crack 
growth  and  arrest  in  a  hybrid  laminate,  the  main  thrust  of  the  investiga¬ 
tion  will  be  to  study  the  behavior  of  the  laminate  as  a  function  of  the 
design  parameters;  fiber/matrix  properties,  buffer  strip  geometry  and 
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deal i zed  Buffer  Strip  Laminate 


the  initial  crack  length.  In  order  to  simplify  the  analysis,  the 
effect  of  any  additional  damage,  e.g.  longitudinal  matrix  yielding 
and  splitting  [9]  or  transverse  matrix  and  fiber  damage  [11]  will  not 
be  considered. 

A  typical  buffer  strip  laminate  usually  contains  angle  plies  as 
well  as  zero  degree  plies.  It  is  felt,  however,  that  much  of  the 
characteristic  behavior  of  the  buffer  strip  region  can  be  represented 
by  the  uni-directional  laminate,  as  a  major  portion  of  the  load  is 
carried  by  these  fibers.  It  appears  that  a  primary  function  of  the 
angle  plies  is  to  prevent  longitudinal  matrix  splitting  in  a  brittle 
matrix  such  as  epoxy.  This  can  be  accounted  for  to  some  degree  in  the 
present  solution  by  allowing  the  matrix  to  support  large  strains  without 
spl itting. 

The  presentation  of  the  solution  will  follow  the  order  of  the 
development  by  the  author  as  this  seems  to  indicate  more  clearly  the 
significant  points  of  the  analysis.  Results  will,  however,  only  be  given 
for  the  final  solutions  corresponding  to  the  buffer  strip  laminate  and 
the  finite  width  strip.  First  the  solution  for  a  crack  in  a  half-plane 
having  arbitrary  shear  stresses  applied  to  the  free-edge,  (which  forms 
the  basis  for  all  the  later  solutions),  will  be  developed.  By  matching 
boundary  conditions  along  the  interface,  the  solution  for  two  different 
adjoined  half-planes  will  be  given  and  then  the  adjoined  half-plane  solu¬ 
tion  will  be  modified  to  account  for  a  second  interface,  resulting  in  the 
buffer  strip  problem.  The  case  of  the  finite  width  strip  with  broken 
fibers  is  obtained  as  a  special  case  of  the  buffer  strip  problem  by  set¬ 
ting  the  shear  stresses  along  the  two  interfaces  to  zero.  These  solutions 
are  all  presented  in  Chapter  II.  Chapter  III  deals  with  the  numerical 
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technique  used  to  solve  a  system  of  linear  algebraic  simultaneous  equa¬ 
tions  coupled  with  a  set  of  linear  Fredholm  integral  equations  of  the 
second  kind.  A  special  method  developed  to  evaluate  accurately  inte¬ 
grals  having  a  cusp  is  emphasized.  Chapter  IV  gives  results  and  con¬ 
clusions  for  the  buffer  strip  problem  and  the  finite  width  strip. 
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CHAPTER  II 


FORMULATION 

Un1 -Directional  Half-Plane  with  Broken  Fibers 

A  uni -directional  array  of  parallel  fibers  with  an  arbitrary  number 
of  broken  fibers  in  the  form  of  a  notch  is  shown  in  Figure  3.  The 
laminate  is  subjected  to  a  prescribed  shear  stress,  T,(y),  along  the 

a 

free  edge  in  addition  to  a  remote  uniform  tensile  stress  in  the  axial 
direction.  Fiber  breaks  occur  along  the  x-axis  (axis  of  symmetry)  and, 
since  the  loading  is  symmetric,  only  the  upper  half  of  the  laminate  is 
considered  in  the  analysis. 

The  fibers  are  taken  to  be  of  much  higher  strength  and  extensional 
stiffness  than  the  matrix  and  therefore  all  of  the  axial  load  is  assumed 
to  be  carried  by  the  fibers  with  the  matrix  transferring  load  by  shear 
stresses  as  given  by  the  classical  shear-lag  assumption.  The  axial  fiber 
stress,  o^(y),  and  matrix  shear  stress,  T^(y),  are  then  given  by  the 
simple  relations 

dv^(y) 

=  Ep  -g^  ,  and 

T^(y)  =  [v^(y)  -  v^_^{y)]  .  (1) 

Where  v^(y)  is  the  axial  displacement  of  the  fiber  n  at  the  loca¬ 
tion  y,  Ep  is  the  Young's  modulus  of  the  fiber,  is  the  equivalent 
matrix  shear  modulus  and  h  is  a  shear  transfer  distance.  Because  of  the 
interference  between  fibers  it  is  unlikely  that  G|^  will  be  the  homo¬ 
geneous  matrix  shear  modulus  or  h  the  actual  fiber  spacing,  and  it  is 
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pointed  out  in  [10]  that  these  values  should  be  determined  experimentally 
for  a  given  laminate.  It  is  also  shown  in  [10]  that  a  single  experi¬ 
ment  giving  the  crack  opening  displacement  as  a  function  of  applied 
load  is  sufficient  to  determine  both  the  equivalent  shear  modulus 
and  the  shear  transfer  distance  h  for  a  particular  laminate,  indepen¬ 
dent  of  the  notch  length.  That  is,  these  parameters  are  material  con¬ 
stants  and  depend  only  on  the  fiber  and  matrix  properties  and  the 
fiber  volume  fraction  but  not  on  the  sizes  of  the  damage. 

By  virtue  of  the  shear-lag  assumption  the  longitudinal  and  trans¬ 
verse  equilibrium  equations  become  decoupled  and  the  fiber  axial  dis¬ 
placements  and  stresses  can  be  obtained  without  solving  the  transverse 
equilibrium  equation.  Therefore,  only  the  equilibrium  equation  in  the 
longitudinal  (axial)  direction  will  be  considered.  With  reference  to 
the  free-body  diagram  of  a  typical  fiber-matrix  region  shown  in  Figure  3, 
the  equilibrium  equations  in  the  longitudinal  direction  is  given  by 

Ap  da  (y) 

1 - ^n+l^y^  ■  ^n^^^  "  °  ’  (2) 

for  all  fibers  except  n  =  0,  and 

Ap  da  (y) 

-f"  — dy —  ^  ^  •  (2) 


Using  the  stress-displacement  relations.  Equation  (1),  in  the  above 
equilibrium  equations,  the  following  set  of  differential -difference  equa¬ 
tions  is  obtained; 


ApEph  d^^ 

,2  Vl 


Got 


dy‘ 


2"n  '^n-l 


0  ,  and 


(4) 


ApErh  d  v„ 

F  F  0  .  /  > 

-  — ^  ''l  "  ''o  "  ^a^^^  • 


dy‘ 


(5) 
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Noting  the  coefficient  of  the  second  derivative  term  in  the  above  equa¬ 
tions,  the  following  changes  in  the  variables  are  suggested: 


let  y 


1 


ApEph 


G»,t 


dv 


"^n  "  dy  ’ 


and 


=  a 
n  o 


1 


Aph 


EpG^t  -n  • 


(6) 


Algebraic  manipulation  then  gives 


a„  =  0 
n  ° 


dV 

n 

-q —  and  =  0 
dn  n  ° 


Oit  ^'^n"'^n-l^  ’ 


(7) 


where,  n  ,  0^  and  V^(n)  are  non-dimensional. 


The  resulting  equilibrium  equations  in  non-dimensional  form  are  given  by 
'f  +  -  2V^  +  =  0  ,  and  (8) 


d2v„ 


dn‘ 


d^V, 


dn 


2""  ''1  ■  "0  “  ’ 


(9) 


where. 


'a<">  ' 


1 


Epth  Tg(y) 


Ar-Glu,  0 
F  M  ° 


These  differential -difference  equations  are  reduced  to  differential 
equations  by  introducing  the  even  valued  transform  as 


V(n,0)  =  x  V^(n)  cos[(n  +  Y)0] 
n=0 


from  which 


TT  - 


V^(n)  =  7  /  V(n,0)cos[(n  +  ^)0]d0 


(10) 


(11) 


Making  use  of  the  above  transformation  and  the  orthogonality  property  of 
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the  circular  functions,  the  two  equilibrium  equations  may  be  written  as 
one  equation  valid  for  all  values  of  n  and  n  as  follows: 

f  /  cos(0)]vj  cos[(n +^)0]d0 

o  7T  -j 

=  7  /  T  (n)cos(0/2)cos[(n +  4-)e]de  , 

IT  0  a  ^ 

which  is  of  the  form 

2  ^  1 

-  /  F(n,e)cos[(n +-2)0]d0  =  0  for  all  n  and  n  . 

0 

Noting  the  definition  of  V(n,0)  in  Equation  (10)  and  (11)  it  is  seen  that 
the  function  F(n,e)  is  even  valued  in  e  and  therefore,  if  the  integral  is 
to  vanish  for  all  the  function  F(n,0)  must  be  zero.  The  single  equa¬ 
tion  specifying  V{r\,Q)  is  then 
2- 

-  6^V  =  ^g(Ti)cos(e/2)  ,  (12) 

dn 

where,  6^  =  2[1-  cos(0)]  4  sin^(0/2)  . 

The  solution  to  the  problem  of  vanishing  stresses  and  displacements 
at  infinity  and  uniform  compression  on  the  ends  of  the  broken  fibers  will 
now  be  sought.  The  complete  solution  is  obtained  by  adding  the  results 
corresponding  to  uniform  axial  stress  and  no  broken  fibers  to  the  follow¬ 
ing  solution.  The  appropriate  boundary  conditions  are  as  follows: 
dV  (n) 

V^(ri)  =  0  ,  — —  =  0  ,  as  n  ^  for  all  fibers,  (13) 

dV^(n) 

— ^  1  5  at  n  =  0  ,  for  all  broken  fibers,  (14) 

V^(n)  =  0  ,  at  n  =  0  ,  for  all  unbroken  fibers.  (15) 

The  complete  solution  to  Equation  (12),  satisfying  vanishing  stresses 
and  displacements  at  infinity,  is  given  by 
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V(n,6)  =  A(6)e-*''  -  J-"  sinh[6(^-t)]^j{t)dt  ,  (16) 

n 

where  the  function  A(9)  is  yet  unknown.  The  remaining  two  boundary  con¬ 
ditions  give 

— ^  /  [-6A(e)  -  cos(0/2)/  cosh(6t)T  (t)dt]cos[(n +i)6]d0  =  - 1 , 

'I''  "  0  0  “  (17) 

for  all  broken  fibers,  and 

V  (0)  =  -  /  [A(e)  +  /°sinh(6t)Tg(t)dt]cos[(n  +  -^)e]d0  =  0  , 

"  0  0 

for  all  unbroken  fibers.  Equation  (18)  is  solved  exactly  by  taking 

f  GO  M  1 

A(e)  +  cos(9/2)_  j.  sinh(6t)^  (t)dt  =  Z  cos[(N*  +  m  +  ^) 0]  ,  (19) 

^0  m=l 

where,  M  is  the  number  of  broken  fibers.  By  eliminating  A(0)  between 
equations  (17)  and  (19),  the  stress  boundary  condition  reduces  to 

O  TT  M  •]  1 

-/  E  B  cos[(N*  +  m  +  y)0]cos[(n  +  ^)0]sd0 

^  0  m=l  ""  ^ 

IT  *1  00  —  <5 1 

+  —  /  cos(0/2)cos[(n  +  i)e]  /  e  t  (t)dt  de  =  1 ,  (20) 

0  ^  0  ^ 

for  n  =  N*  +  1 , . . . ,N. 

For  a  given  shear  stress  distribution,  'rg(t).  Equation  (20)  reduces 
to  a  set  of  linear  algebraic  equations  in  the  Fourier  constants  B^  .  From 
Equations  (16)  and  (19),  A(0)  may  be  eliminated  to  obtain  y{r\,e)  in  terms 
bf  the  constants  B^^  .  Recalling  the  relation  between  V(n,e)  and  V^(n), 
an  expression  can  be  obtained  for  the  axial  fiber  displacement  V|^(ti)  as 

^  -6n  M  11 

V  (n)  =  -  /  e  E  B^  cos[(N*  +  m  +  y)0]cos[(n +-^)0]de 
"  ^  0  m=l  ^  ^ 

_  L  /  f  [){S,r],t)  x,(t)dt  cos[(n +  i)0]d0  ,  (21) 

TT  0  6  0  ®  ^ 

X  -6ln-tl  ~6(Ti+t) 

where  D(6,n5t)  =  e  ^  -  e 
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The  axial  fiber  stress  is  obtained  by  differentiating  Equation  (21)  with 

respect  to  n  and  is  given  by 

dV^(Ti) 

^n  dn 

=  -  -  /  6e  E  cos[{N*  +  m  +  i)0]cos[(n  +  i)0]de 
0  m- 1 

1  IT  «>  f  -i5(n+t)  -6|n-tK  _  1 

-^/cos{0/2)/  |e  -pe  >  Tg(t)dt  cos[(n +^)0]d9, 

(21a) 

where,  p  =  1  for  t  £  n  s  and 
p  =  -  1  for  t  >  n  . 


Adjoined  Uni-Directional  Half-Planes 
Figure  4  shows  two  uni -directional  half-planes  of  different  fiber 
and  matrix  properties  which  are  assembled  to  form  adjoined  half-planes. 
Both  planes  may  have  an  arbitrary  number  of  broken  fibers.  Superscripts 
I  and  II  are  used  to  designate  quantities  corresponding  to  plane  I  and  II 
respectively.  The  normalized  spatial  variables,  n  and  5,  in  the  longi¬ 
tudinal  direction  are  related  by 


/w 


I 


(22) 


The  shear  stresses,  T^(n)  and  tJ^(c),  at  the  interface  are  normalized 

a  a 

with  respect  to  material  properties  of  plane  I  and  II  respectively,  and 
are  related  to  the  actual  shear  stresses  as  follows: 


A 


Epth 

ApGp, 


T^(y) 

a 

r~ 


and 
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Figure 


(23) 


Eptnyl  ^^(y) 

*  oo 

where,  T^(y)  and  T^^(y)  are  the  actual  shear  stresses  at  the  interface 

d  cl 

acting  on  plane  I  and  II.  Since  the  remote  displacements  and  hence  the 
remote  strains  are  assumed  equal,  the  remote  stresses  in  planes  I  and  II 
must  be  related  by 


(5)  - 


00 


00 


(24) 


If  Tj(n)  and  Ti[^(5)  are  known  then  the  solution  for  each  of  the  half- 

a  a 

planes  is  the  same  as  that  developed  in  the  previous  section.  Equation  (.20). 
The  Fourier  constants  and  corresponding  to  broken  fibers  in  plane  I 
and  II  can  then  be  obtained  by  solving  the  following  sets  of  equations 
for  known  Tj(n)  and  tJ^(c): 

a  a 


I'i]  T  o  IT  ]  1 

z  -  /  cos[(N*-,  +m  +  ^)e]cos[(n  +  y)e]6d0 
m=l  0  '  ^  ^ 

oTT  1 

+  ^/  cos(0/2)cos[(n +  ^)9]/  e  Tg(t)dt  de  =  1  , 

^0  0 
Mo 

^  T  T  ?  ^  1  1 

and  E  ^  /  cos[(N*2  +  m +  7)9]cos[(ji  +  ■^)0]6de  , 
m=l  ^  0 

qTT  1  oo”  6txT 

+  -/  cos(0/2)cos[(£  +  ^)0]/  e  ^;(t)dtd0  =  l,  (25) 

'H'  „  ^  „  a 


for  n  =  N*  +  1 ,  +  2, . . .  ,Ni  and  A  =  +  1 ,  +  2 , . . . ,N2  , 

where,  M.j  and  M2  are  the  number  of  broken  fibers  in  planes  I  and  II.  The 
normalized  displacements  of  a  fiber  in  plane  I  and  II  are  then  given  by 

V^(ri)  =  -  /  i:  cos[(N-, -f  m  +  ^)0]cos[(n +  y)0]e  dd 

'•  0  rti=l  I"  I  ^ 

-  -  /  cos[(n+l)0]  r  0(6, n,t)  ^ht)dt  de  ,  (26) 

0  ^  “^0  ® 
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When  the  above  two  half-planes  are  joined  together  the  shear  stress 
along  the  interface  is  unknown,  but  from  equilibrium  the  shear  stresses 
T^(y)  and  T^^(y)  acting  on  each  of  these  two  half-planes  must  be  equal 
and  opposite.  Further,  as  the  shear  stress  is  directly  related  to  the 
distortion  of  the  matrix  from  the  shear-lag  assumption,  it  follows  that 
these  stresses  must  be  proportional  to  the  difference  in  the  displacement 
of  the  first  fibers  of  plane  I  and  II.  These  conditions  result  in  the 
following  two  equations: 

Tgly)  =  -  'rg^(y)  5  and  (29) 


TgCy)  =  (G^/h)' 


(30) 


where,  (Gj^/h)^  is  the  equivalent  shear  stiffness  of  the  interface.  It  is 
interesting  to  note  that  in  a  continuous  elasticity  solution  one  would 
match  surface  tractions  and  displacements  at  the  interface  while  in  the 
present  discrete  modeling  solution  the  shear  stresses  are  required  to  be 
equal  and  the  shear-lag  relation.  Equation  (30),  takes  the  place  of  the 
displacement  equality. 
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Substituting  for  the  actual  displacements  in  terms  of  normalized  values 
using  Equation  (28),  and  recalling  the  relation  between  the  actual  and 
normalized  shear  stresses  from  Equation  (23),  Equations  (29)  and  (30)  re¬ 
duce  to 


^3^(5)  =  -  (G|v|/h)^(h/G|^)^^ 


T3(n)  =  vj(n)  -  Ri  vj^s)  (G^/h)V(G^/h) 


where . 


ApEph 


ApEph 


Using  Equations  (26)  and  (27),  vj(n)  and  vj^(6)  can  be  found  and 
substituted  into  Equation  (32)  thus  resulting  in  an  integral  equation 
for  ^g(n)  in  terms  of  the  Fourier  constants,  and  ,  and  the  nor¬ 
malized  shear  stresses,  Tg(n)  and  '^^^(S)-  The  spatial  variable  6  and 
the  normalized  shear  stress  tJ^(0  may  be  eliminated  from  the  above  in- 

a 

tegral  equation  and  also  from  the  equations  for  the  Fourier  constants. 
Equation  (25),  resulting  in  the  following  set  of  governing  equations; 
two  series  equations  and  one  linear  integral  equation: 


E  B^  -  /  cos[(Nt  +  m  +  i)e]cos[(n+i)0]6de 
m  TT  I  ^  ^ 


m=l  ^  0 


O  TT  1  ”  -T 

+  -  /  cos(e/2)cos[(n +  y)e]  /  e  t  (t)dt  de  = 


E  B^  -  /  cos[(N2  +  m  +  ^)0]cos[(£  +  -2)e]5de 
m=l  m  ^0 


"^1^  I, 


cos(0/2)cos[(£  +  i)0]  f  e  T  (t)dt  de  - 


^  0 
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for  n  =  +  1 .  ,Ni  and  +  1 » •  •  •  5N2  ,  and 

TT  T  I 

T^(n)  =  G.n  -/  s  cos[(Nt +  ni  +  i)0]cos(e/2)e  de 
a  0  m=l  ‘  ^ 


K 


-  R  G 

I  I  I  IT 


2  II  .1  -5in 


il  f  2  £  cos[(N*  +  in  +  l)«> 


de 


0  m=l 


0  ”  0 


2  /  cos^( e/2) 
6 


D(5,n,t) 


where,  G^i  =  (Gj^/h)^  (h/G^^)  ^  ,  G12  -  (G^/h)  ^(h/Gj^)^^  and  61  -  S/R-j 


de  T:(t)dt 

a 


(36) 


The  Buffer  Strip  Laminate 

Figure  5  shows  a  finite  width  strip  (buffer  strip)  between  two  half¬ 
planes  of  different  or  of  the  same  fiber  and  matrix  properties.  One  of 
the  half-planes,  plane  I,  and  the  finite  strip,  region  II,  may  contain  an 
arbitrary  number  of  broken  fibers.  The  normalized  spatial  variables  tij  5 
and  c  are  related  to  each  other  in  the  same  way  as  in  Equation  (22).  The 
normalized  shear  stresses  Tg(n),  t^^(5)  and  t^^^(^)  are  related 

to  their  corresponding  actual  shear  stresses  in  the  same  manner  as  in 
Equation  (23).  Further,  the  remote  strains  in  all  three  regions  are 
assumed  equal. 

The  solutions  for  planes  I  and  III  can  be  obtained  for  known  T^(ri) 
and  the  half-plane  solution  and  therefore  we  need  to  de¬ 

termine  the  solution  for  the  finite  strip,  region  II,  only.  This  solu¬ 
tion  is  developed  by  considering  the  half-plane  shown  in  Figure  3  with 
the  following  special  condition.  It  may  be  assumed  that  the  matrix  be¬ 
tween  the  fiber  Nw  and  (Nw  +  1)  splits  all  the  way  to  infinity  and  that 
an  external  shear  stress  T^^(y)  is  applied  on  the  surface  of  the  split 
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as  shown  in  Figure  6.  The  introduction  of  a  split  with  or  without  the 
external  applied  shear  stress  results  in  two  special  equilibrium  equa¬ 
tions  for  the  two  fibers  on  either  side  of  the  split  [9,  10,  11]  in 
addition  to  those  valid  for  the  free- edge  fiber  and  the  generic  fiber, 
n,  given  by  Equations  (8)  and  (9).  These  equilibrium  equations  in  the 
normalized  form  are  as  follows: 


H2\/n 

d^2  1  0  a 


for  fiber  0  , 


d2yll 

*  Cl  -  K "  '"i  =  0  .  f"--  "  . 

dS 


.^2  ^Nw  Nw-1 

dS 


-  ,  for  fiber  Nw  ,  and 


d2vll 

"^2  -  Cl  ' 

The  left  hand  sides  of  Equation  (39)  and  (40)  can  be  reduced  to  the 
standard  form  by  adding  or  subtracting  a  term  on  both  sides. 

Making  use  of  the  transforms  similar  to  Equation  (10)  and  (11)  and  fol¬ 
lowing  the  procedure  of  the  half-plane  problem,  the  single  equation 
specifying  the  transformed  normalized  displacement  V^^(5,e)  can  be 

obtained  as 
2"  1 1 

dj/  =  ^J^(c)cos(0/2)  +  [g(c)  -  ri^(?)]F^  ,  (41) 


where, 

g(5)  =  ,  and  F^  =  cos[(Nw  +  ^)0]  -  cos[(Nw  +  |)0j  . 

The  solution  to  Equation  (41)  satisfying  vanishing  stresses  and  dis¬ 
placements  at  infinity  and  unit  compression  on  the  crack  surface  is  then 
gi  ven  by : 
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applied  shear  stress 


M2 

Z  ^1}  —  f  cos[(N^  +  m  +  i)0]cos[(n  +  i)0]d0 
ni=  1  0 

9  ^  1  “  ”<5t  _JT 

+  -/  cos(e/2)cos[(n +  ^)0]  /  e  T“(t)dt  d0 

+  — /  cos[(n+i)0]  /  e  fg(s)  -  ds  do  =  1,  (42) 

^0  ^  0  I  J 

for  all  broken  fibers. 

If  T^^(s),  g(t)  and  tJ^(s)  are  known.  Equation  (42)  reduces  to  a  set  of 

d  D 

linear  algebraic  equations  in  and  can  be  solved  directly. 

When  the  above  finite  strip,  region  II,  is  introduced  between  the 

two  half-planes,  planes  I  and  III,  as  shown  in  Figure  5,  it  results  in 

two  interfaces  each  similar  to  the  one  discussed  in  the  problem  of  adjoined 

half-planes.  The  interfacial  shear  stresses  ■t^(n),  t^^(^),  Tj^(g)  and 

d  a  D 

tP^(?)  can  be  obtained  in  the  same  manner  as  in  the  previous  solution, 
b 

that  is,  using  the  following  relations: 

Tg(y)  =  -  Tg^(y),  and 

Tg(y)  =  (G,v,/h)''^  [vj(y)  -  vj^(y)],  between  planes  I  and  II,  (43) 
and 

T^^(y)  =  T^^^(y),  and 

'rJ^(y)  ="(GM/h)^^  [v^^(y)  -  v^^^(y)]  between  planes  II  and  111,(44) 
b  Nw  0 

where  (G^^/h)^^  and  (Gj^/h)^*^  are  the  shear  stiffnesses  of  the  respective 
interfaces.  By  definition  g(§)  is  given  by 

"  '^Nw+l  ■  '^Nw  • 

Equations  (43)  -  (45)  along  with  the  two  equations  for  the  Fourier  con¬ 
stants  from  the  stress  boundary  conditions  on  the  broken  fibers  in  region 
I  and  II,  Equations  (34)  and  (42),  represent  the  complete  solution  for  the 
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finite  width  buffer  strip  problem.  In  Equation  (45)  it  must  be  noted 
that  is  the  normalized  displacement  of  fiber  Nw  of  the  finite  strip, 
whereas  is  the  normalized  displacement  of  the  first  fiber  of  a  new 

half-plane  having  applied  shear  stresses  along  the  free  edge  and 

fiber/matrix  properties  of  the  finite  strip. 

As  in  the  previous  problem,  first  the  actual  displacements  are  sub¬ 
stituted  in  terms  of  their  corresponding  normalized  displacements,  then 
the  normalized  displacements  are  substituted  in  terms  of  the  normalized 
shear  stresses  and  the  Fourier  constants.  Further,  and 

can  be  eliminated  knowing  the  relationship  between  the  actual  shear 
stresses.  Equations  (43)  and  (44).  Then  the  final  set  of  governing  equa¬ 
tions  for  the  finite  width  buffer  strip  problem  can  be  listed  as  follows; 

Ml 

I  T  Q  TT  1  1 

E  B;:  -  /  cos[(Nt  +  m  +  4-)9]cos[(n  +  4-)e]6de 

in=l  "  '  0  '  2'  ^  2'  J 

pTT  1  co-6t^j 

+  — /  cos(0/2)cos[(n +4-)e]  f  e  T^(t)dt  de  =  1  ,  (46) 

^0  ^  0  ^ 


E  —  /  cos[(Ni  +  m  +  i) 0]cos[(£  +  i) 0]6de 


^17  9^  1  “~^ll'_T 

-n —  /  cos(e/2)cos[(£  +  ^)0]  /  e  TVt)dt  d0 

•^4-  7T  c  a 


+  I  /  F^cos[(£  +  i)e]  /  e  fg(s)  -  ^J^(s)]  ds  do  =  1 


for  n  =  N-|  +  1 ,. . .  ,Ni  and  £  =  N2  +  1 , . . .  ,N2  , 

_T  ,  '  T  pi:  .  n  “on 

Tg(n)  =  E^  G^.-j  cos[(Nlj'  + m  +  '2-)0]cos(0/2)e  do 


^  TT  ?  *  1  “'^in 

^  ~ cos[(N5  +  m  +  y)0]cos(0/2)e  do 

^^1  m  1  I  I  ir  0  ^  ^ 
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oo 


and 


G  TT  2 

^  /  cos  W2).  [0(6, n,t)  +  G^2  0(5-,  ,n,t)]de  ;^(t)dt 


-  / 

0  0 

+  /“g.-,  —  /  cos_(|/2j_p2  0(6, n/R^  ,s)d0  [g(s)  -  ^J^(s)]dS: 


0  .  ^  0 


(48) 


'"'2  -65 

g(c)  =  -  ^  -  /  cos[(N^-i-ni-f-i)e]cos[(Nw-i-i)0]e  de 

m=l  0 


,  ^12  1  /  cos(e/2) 

~  J  n  ^  J 


^  cos[(Nw-i-|-)0]  D(6,t/Rp5)d0  T^(t)dt 

1 


-I, 


0  R^^  0 


CO  1  TT  COS[(Nw 9 
-I-  /  -  /  - j — ^ - F  D(6,5,s)d0  g(s)ds 


b 


0  0 


00-1  TT  I  ___ 

J.  ]_  J.  J  COS 
0  0 


2(,/2)  F^cos[(Nw4)9]'] 


D(6,5,s)de  T^^(s)ds, 


(49) 


M 


2  ”  u  4, 

^  Gn9  ~  T  cos[(N^ -l■nl■^i)0]cos[(Nw-l-i)0]e  d0 

1^1  /  coi(^  cos[(Nw  +  i)9]D(6,t/Rp5)d0  T3(t)dt 
R-j”^  0 

G.2  Tr'cos[(Nw  +  ^)0]  2 


•65 


/ 


CO  Gj  oG 

-  / 

0  R 


F  D(6,5,s)d9  g(s)ds 


where 


00  G-i  T7 


'"23  cos^(0/2) 


0  L'2 


0(623533) 


cos[(Nw-t-i)0]  o 
- -  F  0(6,5, s) 


d0  Tj^(s)ds  (50) 
b 
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-  S/R-j  5 

62  “  ^ 

^11  =  (Gjy|/h)^'^(h/Gn^)^  5 

Gi2  =  (G^/h)''^(h/G^)^^  ,  and 

G23  =  (G^/h)^^(h/Gf^)^^^  . 

The  solution  reduces  to  two  series  equations  and  three  linear  integral 
equations . 


Uni-Directional  Finite  Width  Strip  with  Broken  Fibers 

In  this  section  a  solution  is  obtained  for  a  uni-directional  finite 

width  strip  containing  an  arbitrary  number  of  broken  fibers  as  a  special 

case  of  the  solutions  already  developed.  With  reference  to  Figure  5  it 

II  II 

can  be  seen  that  if  the  interfacial  shear  stresses  (y)  and  (y) 
acting  on  the  region  II  are  zero  then  it  results  in  a  finite  width  strip 
with  broken  fibers  subjected  to  a  uniform  remote  axial  stress  as  shown 
in  Figure  7.  Therefore  the  solution  to  the  finite  width  strip  problem 
can  be  obtained  as  a  special  case  of  the  solution  developed  for  region 
II  of  the  buffer  strip  laminate  by  setting  and  ^5^(2)  to  zero  in 

Equations  (42).  With  reference  to  Figure  7  these  equations  are  given  as 


I'l  o  TT  +1  1 

E  ^  /  cos[(N  +  m  +  ^)e]cos[(n +  ^)e]d0 

«  TT  O  -I  ” 

+  ^  /  F  cos[(n+^)0]  /  e  g(t)dt  de  =  1  , 

0  0 

for  all  broken  fibers,  and 

M  0.1  l-SC 

g(5)  =  -  E  ^  ^  cos[(N*  +  m  +  ^)0]cos[(Nw  +  ^)9]  e  de 
m  =  1  ^0 

CO  ^  IT  cos[(Nw  +  y)e]  „ 

+  /  -^  /  - g— ^ - F^  D(6,c,t)d9  g(t)dt 


(51) 


(52) 
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where 


=  cos[(Nw  +  ■^)e]  -  cos[{Nw  +  |)e]  ,  and 
"  '^Nw+1  "  ^nvi  ° 

Since  by  definition  is  the  displacement  of  the  first  fiber  of 

a  new  half-plane  having  fiber/matrix  properties  of  the  finite  width 
strip  and  no  applied  shear  stress  along  the  free  edge^  it  must  be  equal 
to  zero.  Therefore  g(5)  is  given  by 

g(5)  =  -  Vu„  .  (53 

Therefore  the  solution  of  the  finite  width  strip  reduces  to  one 
series  equation  coupled  with  one  linear  integral  equation.  The  location 
and  the  number  of  broken  fibers  are  arbitrary  except  that  the  fiber 
breaks  must  be  along  the  x-axis.  The  solution  obtained  holds  for  a 
central  notch,  an  edge  notch,  an  off-center  notch,  or  for  multiple 
notches  along  the  x-axis. 


28 


CHAPTER  III 


SOLUTION  TECHNIQUE 

In  all  the  problems  dealt  with  so  far,  the  solution  reduces  to  one 
or  two  series  equations  coupled  with  one  or  more  linear  Fredholm  in¬ 
tegral  equations  of  the  second  kind.  Since  there  is  no  exact  closed 
form  solution  available  to  solve  such  a  system  of  equations,  a  numerical 
procedure  is  developed  and  presented  in  this  section.  The  technique 
makes  use  of  a  method  given  by  Riez  [12],  to  solve  a  linear  Fredholm 
integral  equation  of  the  second  kind  defined  within  a  semi -infinite 
interval  of  integration.  The  solution  is  based  on  the  fact  that  a  given 
integral  over  a  semi-infinite  interval  may  be  approximated  by  the 
Gauss-Laguerre  quadrature  rule  as 

CO  k  -Xn 


/  f(x)dx  =  z  w.  e  f(x^ ) 
0  i=l 


(54) 


where  x^-  is  the  i^*^  zero  of  the  Laguerre  polynomial,  L|^(x) ,  and  w^-  is 


the  corresponding  weight  function  given  by 
X.I 

w.  = 


1 


^  [(k+1)  L^^^(x.)]" 

The  Laguerre  polynomial  L|^(x)  is  given  by 

I  ,  ,  X  d'^(x^e"^) 

L.  (x)  =  e  — ^ . 

dx^ 


(55) 


(56) 


Since  the  form  of  the  equations  for  each  of  the  solutions  is  the 
same,  the  development  and  application  of  the  numerical  procedure  can  be 
demonstrated,  without  loss  of  generality,  by  taking  the  equations  corre¬ 
sponding  to  one  of  the  solutions.  Consider  then  the  solution  corre¬ 
sponding  to  the  problem  of  adjoined  half-planes  given  by  Equations 


I  CO  -I  1 

Z  -  /  cos[(Nt  +  m  +  i)e]cos[(n  +  -p)e]6d0 

m=l  ^  0  '  ^  ^ 

+  -  /  cos(0/2)cos[(n +  l)e]  /  e  Tht)dt  de  =  1  ,  (57) 

'n"  ^  ^  n  ^ 


M2 

z  -  7cos[(N5  +  m-i-i)e]cos[{j’ +I)e]6d0 
T  m  TT  ^  ^  ^  ^ 

m=1  0 

-  — /  cos(e/2)cos[(£  +  i)0]  /  e  T^(t)dt  de  =  1  ,  (58) 

’  0  “^0 

for  n  =  +  1 5  a « 0  jN-j  and  £  =  N2  1 9 •  5^2  9  snd 


1 

x(n)=G.-,-/  E  B^  cos[(n:,' +in  +  ^)e]cos(0/2)e  de 
^  ^  0  m=l  ' 

M2 

-R  G  i  /  E  B^^  cos[(N,  +  m+ i)e]Gos(e/2)e 
I  '  I  tt  1  rn  c  c, 

0  m=l 


2 

D(6,n,t)  +  D(5pn,t) 

X  de  7(t)dt  .  (59) 
a 

Integrals  over  the  spatial  variable,  t,  in  the  above  equations, 
defined  over  the  semi -infinite  interval,  can  now  be  replaced  by  the  above 
series  representation  to  yield 


r  1  r'^  COS  (6/2) 

\3  ^  -f  J  J  o 

0  0  ^ 


|Vll 

E  { 
m=l 


-  /  cos[(N*  +  m  +  ^)0]cos[(n  +  ^)0]6de  |  bJ, 
^  0 


"2 
E  < 
m=l 


klpiT  1  “'^^i  ^i 

+  E  i  -  /  cos(e/2)cos[(n +'^)0]e  de  w^e 
j=1  0 

—  /  cos[(N2  +  m  +  -^)0]cos[(£  +  ^)e]6d0  i  B^ 


(60) 
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-  I.  \ -  /  cos(e/2)cos[(£  +  i)e]e  ^  ^  do  w-e  - 1 , 

i  (61 

for  n  =  N*  +  1 , . . .  ,Ni  ,  and  £  =  +  1 , • • • >N2  ,  and 


ij-  K(ni,t.)w. 


M 


-6n. 


m=l 

M2 
+  i:  ( 
ni=l 


E  G-i  I  /  cos[(N^  +  ni  +  ^)e]cos(e/2)e  ^ 


de 


RlGi^  ^  /  cos[{N2  +  m  +  ^)0]cos(0/2)e  ‘  '  de 

i  =  1 .  ,k 


X  b“  =  0. 


(62) 


where,  6^^ 


=  1  for  i  =  j 
=  0  for  i  5^  j  , 


and  K(Ti.,tj) 


^*^'1  /  cos^(e/2) 


'12 


[D(6,n^,tj)  + -^  D(6i,n^-,tj)]de  . 


Therefore,  the  solution  of  the  integral  equation  coupled  with  a  set  of 
series  equations  reduces  to  solving  a  system  of  linear  algebraic  equations 
in  and  explicit  values  of  T^(n)  at  the  quadrature  points.  The 

above  system  of  equations  is  solved  by  the  method  of  Gauss-elimination 
with  partial  pivoting. 

It  must  be  pointed  out  that  the  terms  in  the  integral  equation  of 
the  form  have  cusp  at  n  =  t  and  the  quadrature  representation  in 

the  standard  form.  Equation  (54),  results  in  large  inaccuracies  especially 
for  large  values  of  n*  This  has  been  overcome  by  modifying  the  above 
quadrature  rule  by  deleting  selected  terms  in  the  series  and  replacing 
them  by  closed  form  integration  in  the  neighborhood  of  the  cusp.  With 
reference  to  Figure  8,  the  integral  over  the  semi -infinite  interval  is 
then  represented  by  the  modified  equation 
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k  f 

/  K(n.5t)  Tj(t)dt  =  s  K(n.5t.)  xHtOe  ^  w. 
0  ^  3  j=]  1  J  a  J  J 


t,- 


+  K(n^,t.^-,)  Tg(t.^l)e 


^i+1 


where. 


ti 


-I, 


ti^-|+A2 


-I/. 


/  K(n^,t)  Tg^(t)dt+  /  K(n^,t)  Tg(t)dt 

‘  t  • 


(63) 


=  (t^.i  -  ■t^._2)/2  ,  and 

^2  ""  ^^i+2  "  ' 


Since  Xg(t)  is  yet  unknown,  the  two  integrals  on  the  right  hand  side 
of  Equation  (63)  can  not  be  evaluated  in  a  closed  form.  However,  if  the 
function  T^(t)  is  taken  to  be  the  average  of  the  terminal  values  within 
each  of  the  intervals,  the  above  two  integrals  can  be  expressed  as 


f  KCn^^t)  g(t)dt  =  2 - /  K(n.i,t)dt  , 

S'-l'^l  S'-r^i 


and 


ti^l+A2 


f  K(rn,t)  T^(t)dt  = 


-I .  V‘2 


‘i 


1’  '  a' 


/  K(ni,t)dt  .  (64) 

L  I 


'1 


Substituting  Equation  (64)  into  Equation  (63)  and  rewriting 
results  in 


“  T  ^  T  ^i 

/  K(n.,t)  Tj(t)dt  =  z*  K(ni,t)  ^^(tje  ^ 

0  '  ^  j=l  1  a  -1 

t 


J ' "  j 


1  ^  T 

+  2  f  K(n^,t)dt  x^(t._.|) 


^i-r^l 
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t.^l+A^ 


/  K(n^-,t)dt  +  /  K{Ti^-5t)dt  ' 


1  ^i+1'^^2  j 

+  2  f  K(Ti^,t)dt  Tg(t^-)  , 


where  E  excludes  terms  corresponding  to  j  =  i  -  1 ,  i  and  i  +  1 . 

The  quadrature-rule  in  the  modified  form.  Equation  (65),  is  used  to 
represent  those  integrals  in  which  the  integrand  has  terms  with  cusp. 

The  two  series  equations  remain  unaltered  and  the  substitution  of  Equa¬ 
tion  (65)  into  Equation  (59)  results  in 


k  t 

E*K(„i,t.)w.  e  ^  +  a._,  T-J(t,_,)  +  (Hapj'(t.)+a.^, 

J  * 

Ml  C  ^ 

oTT  ^  1  T 

-  z  “/  cos[(Ni  +  m +  -^)9]cos( 9/2)e  d9  ^ 

m=l  '  ^  0 


“2  2  ^  ,  -<S-,n^- 

+  E  |R-|G^.1  —  /  cos[(N5  +  m  +  ^)9]cos(e/2)e  do 


X  =  0 
m 


i  =  1 , . . .  jk 


where. 


a.  ,  =  I  /  K(ni,t)dt 

"  t  -A 

^i-1  ^1 


t.^l+A2 


1  ■■  1  ^ 

•  = -i-  /  K(n4,t)dt  +  y  /  K(n,-,t)dt  ,  and 

I  ^  4.  A  •  4>  ' 


ti-r^i 


,  ^i+l'^^Z 

a^+l  =  ^  /  K(n^.at)dt 

^i 


This  again  reduces  to  solving  set  of  linear  algebraic  equations  as  before. 
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CHAPTER  IV 

RESULTS  AND  CONCLUSIONS 

First,  typical  results  are  given  for  a  Graphite/ epoxy  laminate  con¬ 
taining  a  buffer  strip  as  shown  in  Figure  9.  The  crack  growth  behavior 
of  the  lamina  is  studied  by  varying  the  buffer  strip  material,  the  width 
and  the  thickness  of  the  strip  and  the  initial  crack  length.  The  three 
materials  considered  for  the  buffer  strip  are  Nylon,  Kevlar  and  S-glass, 
all  in  an  epoxy  matrix.  The  matrix  and  fibers  are  assumed  to  be  linearly 
elastic.  The  failure  criterion  is  that  a  fiber  fails  upon  reaching  its 
ultimate  failure  stress  as  determined  from  an  unnotched  coupon  test. 

Figure  10  presents  results  corresponding  to  initial  crack  growth  in 
plane  I,  crack  arrest  at  the  interface,  crack  growth  in  the  buffer  strip 
and  subsequent  lamina  failure.  In  these  results  all  fibers  are  of  the 
same  cross-sectional  area  and  in  all  cases  the  buffer  strip  is  ten 
fibers  wide.  Since  all  the  buffer  strip  candidate  materials  are  of 
lower  modulus  than  that  of  the  parent  laminate,  which  in  this  case  is 
Graphite/epoxy,  the  stress  concentration  factor  at  the  near  end  of  the 
notch  (nearer  to  the  interface)  is  always  higher  than  that  at  the  far 
end  of  the  notch.  The  solid  line  in  Figure  10  represents  the  remote 
stress  required  to  initiate  crack  extension,  (fail  the  first  unbroken 
fiber  in  front  of  the  notch,  fiber  A  in  Figure  9).  The  remote  stress 
required  to  fail  the  lamina  catastrophically,  (fail  the  first  fiber  in 
plane  III,  fiber  B  in  Figure  9)  is  given  by  the  broken  line  in  Figure  10. 
Both  these  stresses  are  functions  of  the  initial  crack  length  and  decrease 
with  increasing  length.  The  crack  growth  takes  place  by  breaking 
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Figure  9.  Buffer  Strip  Laminate  with  Initial  Damage 


consecutive  fibers  from  the  crack  tip  to  the  interface.  Then,  depending  on 
the  stress  level  required  to  run  the  crack  to  the  interface  and  depending 
on  the  buffer  strip  material ,  the  crack  may  arrest.  It  is  very  interesting 
to  note  that  all  three  buffer  strip  materials  require  an  increasing  stress 
to  continue  the  crack  growth  in  the  buffer  strip,  although  Kevlar  will 
only  arrest  a  crack  if  it  initiates  under  fairly  low  load,  i.e.  initially 
close  to  the  interface. 

Total  lamina  failure  will  occur  when  either  the  buffer  strip  is  fully 
broken  or  when  the  first  fiber  in  plane  III,  fiber  B,  fails.  In  both 
cases  continued  crack  growth  is  unstable  once  fiber  B  fails.  For  the 
particular  lamina  of  Figure  9,  all  fibers  in  the  Kevlar  buffer  strip  fail 
before  fiber  B  attains  its  failure  stress,  whereas  for  S-glass  and  Nylon, 
fiber  B  fails  when  there  are  still  some  fibers  left  unbroken,  i.e.,  the 
crack  jumps  the  buffer  strip.  The  ultimate  lamina  failure  stresses  are 

a  =  0.272  a  for  Kevlar 

“  ul  t 

=  0.395  for  S-glass,  and 

=  0.444  for  Nylon, 

where  ‘is  the  ultimate  fiber  stress  as  determined  from  an  unnotched 

Graphite/epoxy  laminate.  The  material  properties  used  for  these  results 
are  given  in  Table  1.  The  results  of  Figure  10  indicate  Nylon  to  be  the 
best  of  the  three  materials  but  this  is  only  true  if  the  matrix  can 
support  the  very  large  failure  strains  (about  20%)  of  Nylon.  In  a  typical 
angle-ply  laminate  with  Nylon  buffer  strips  and  with  continuous  +45 
graphite  plies,  high  strain  levels  certainly  cannot  be  reached  before 
failing  the  angle-plies  and  continuing  the  crack.  This  behavior  was 
observed  by  Verrette  and  Labor  [4].  The  extension  of  this  study  to 
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account  for  angle  plies  is  surely  necessary  to  represent  accurately  such 
low  modulus  fibers. 


TABLE  1.  FIBER  PROPERTIES 


fiber 

mass  density 
Gms/m^ 

modulus 

MPa 

ult.  stress 
MPa 

Nylon 

84 

2000 

233 

Kevlar 

106 

111400 

2020 

S-glass 

199 

101000 

2800 

Graphite 

155 

300000 

2800 

In  Figures  11-13  the  effect  of  buffer  strip  width  on  crack  growth 
through  the  strip  is  indicated.  The  ultimate  failure  stress  of  the 
lamina  as  a  function  of  buffer  strip  width  is  plotted  in  Figure  14. 

From  Figure  14  it  is  seen  that  for  Nylon  the  ultimate  failure  stress 
continues  to  increase  with  an  increase  in  width,  whereas  for  S-glass  and 
Kevlar  about  ten  fiber  widths  is  optimum. 

Table  2  gives  a  comparison  of  the  ultimate  lamina  failure  stress 
for  a  buffer  strip  of  ten  fibers,  first  for  equal  area  and  then  for 
equal  weights  of  the  fibers.  These  results  are  normalized  with  respect 
to  S-glass.  A  Kevlar  region  of  two  plies  is  approximately  equivalent  to 
one  ply  of  S-glass  in  weight  and  it  is  seen  that  S-glass  is  still  superior 
to  Kevlar.  This  agrees  with  results  obtained  by  Poe  and  Kennedy  [5]. 

Next,  results  are  presented  for  the  finite  width  strip  with  broken 
fibers  forming  a  central  notch  and  subjected  to  a  uniform  remote  axial 
stress  as  shown  in  Figure  7.  The  stress  concentration  factors  for  various 
notch  lengths  (number  of  broken  fibers,  M)  are  obtained  corresponding  to 
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TABLE  2.  COMPARISON  OF  LAMINA  FAILURE  STRESS 


f  i  ber 

equal  area 

^~'^°ul  t 
equal  weight 

Nylon 

0.444 

0.472 

Kevl ar 

0.272 

0.318 

S-gl ass 

0.395 

0.395 

two  widths,  Nw  =  24  and  48.  Following  Hedgepeth  [6]  the  stress  concentra¬ 
tion  factors  for  a  uni -directional  infinite  region  are  obtained  using  the 
relation 

'^N+1  _  4-6-8-  .....  .(2M+2) 

"^N+l  a  3-5-7 . (2M+1) 

where 

N+1  =  index  of  the  first  unbroken  fiber  at  the  notch  tip, 

M  =  total  number  of  broken  fibers  in  the  notch, 

=  stress  concentration  in  fiber  N+1, 
a|^^.|  =  axial  stress  in  fiber  N+1,  and 
=  applied  remote  stress. 

By  comparing  the  above  two  stress  concentration  factors  corresponding 
to  a  given  notch  length  the  finite  width  correction  factor  for  a  uni¬ 
directional  strip  is  obtained  and  is  given  in  Table  3.  Also  given  in 
Table  3  are  the  corresponding  finite  width  correction  factors  for  an 
isotropic  strip  obtained  from  the  following  relation  [13] 

Y2  =  1.0  +  0.1282  (2a/W)  -  0.2881  (2a/W)^  +  1 .5254  (2a/W)^,  (68) 

where, 

Y2  =  finite  width  correction  factor, 
a  =  half-crack  length,  and 
W  =  width  of  the  strip. 
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From  Table  3  it  is  seen  that  the  finite  width  correction  factors 
for  the  isotropic  strip  and  the  uni-directional  strip  are  not  signifi¬ 
cantly  different  for  aspect  ratios  (notch  length  to  strip  width  ratios) 
less  than  0.5.  For  higher  values  of  aspect  ratios  there  is  a  considerable 
difference  between  the  two  cases.  An  aspect  ratio  of  0.9167  corresponds 
to  22  broken  fibers  in  the  24  fiber  strip  and  all  the  applied  load  is 
carried  by  the  two  remaining  intact  fibers  giving  a  stress  concentration 
factor  of  11.0,  (22x1.0/2).  Similarly,  the  other  limiting  case  corre¬ 
sponding  to  the  aspect  ratio  of  0.9583  in  which  46  fibers  are  broken  out 
of  48,  leaving  the  two  end  fibers  to  carry  all  the  applied  load  giving  a 
stress  concentration  of  23.0  (46/2).  These  two  limiting  cases  are  pre¬ 
dicted  accurately  by  the  above  solution  indicating  the  validity  of  the 
model . 

TABLE  3.  FINITE  WIDTH  CORRECTION  FACTORS 


j.^ —  2a  — v| 


2a  =  M  X  d 
W  =  Nw  X  d 
d=fiber  spacing 


2a/W 

Fini te 

width  correction  factors 

uni -directional  strip 

Tcntrnnic 

Nw  =  24 

Nw  =  48 

0.1667 

1 .01334 

1.01125 

1 .02042 

0.25 

1 .02890 

1.02514 

1.03788 

0.3333 

1.05207 

1.04581 

1.06722 

0.5 

1.13040 

1 .11598 

1 .18275 

0.75 

1 .44760 

1 .39280 

1.57762 

0.9167 

2.57412 

3.7 

2.05049 

0.9583 

3.77560 

2.20070 
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From  tho  abovG  rGsults  and  discussions  the  following  conclusions 
can  be  made.  The  method  predicts  the  fracture  behavior  of  a  hybrid 
laminate  in  terms  of  material  properties ^  geometry  and  initial  crack 
length.  The  results  agree  well  with  those  obtained  experimentally  for 
buffer  strips  of  high  modulus  fibers  such  as  S-glass  and  Kevlar,  where 
the  stiffness  of  the  angle-plies  is  very  small  compared  to  that  of  zero 
degree  ply.  For  low  modulus  (high  failure  strain)  buffer  strip  materials 
such  as  Nylon,  in  which  the  stiffness  of  continuous  angle-plies  is  com¬ 
parable  to  that  of  the  axial  buffer  strip  fibers,  the  model  is  inadequate 
due  to  the  basic  assumption  made  in  idealizing  the  laminate  as  a  uni¬ 
directional  composite.  The  method  predicts  the  best  buffer  strip  material 
to  be  one  with  a  low  modulus  and  as  large  an  ultimate  strength  as 
possible,  e.g.,  S-glass  is  superior  to  Kevlar. 

As  stated  earlier,  the  main  aim  of  this  work  was  to  understand  the 
basic  mechanism  of  crack  growth  and  arrest  in  hybrid  laminates,  keeping 
the  model  as  simple  as  possible.  This  goal  has  been  achieved.  However, 
in  order  to  represent  an  actual  buffer  strip  panel  more  realistically, 
the  model  certainly  needs  and  has  the  potential  for,  further  modifications. 
An  immediate  extension  is  to  model  the  panel  as  a  periodic  uni -directional 
hybrid  laminate  simulating  the  regular  placement  of  the  buffer  strips. 

Next,  the  role  of  angle  plies  must  be  accounted  for  by  either  adding  a  cover 
sheet  over  the  uni-directional  laminate  or  by  some  other  means.  Finally, 
the  effect  of  matrix  damage  in  the  form  of  longitudinal  yielding  and 
splitting  at  the  crack  tip  and  at  the  interfaces  must  be  incorporated. 
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